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Abstract. By imposing natural geometrical and kinematical conditions on a conformal 
Killing vector in Bianchi I spacetime, we show that a class of axisymmetric metrics admits a 
O ■ conformal motion. This class contains new exact solutions of Einstein's equations, including 
^ ■ anisotropic radiation universes that isotropise at late times. 
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1. Introduction 

A conformal Killing vector (ckv) £ M generates a local G\ of conformal motions of the metric: 

ds 2 = g^dx^dx" — > e*c?s 2 <^ 

= Cd a g^u + 9audf,C + g^duC = 2ipg^ (1) 

where ^ is the conformal factor, with ?/> = if £ M is a Killing vector (kv). Conformally 
flat spacetimes admit a maximal Gi$ of conformal motions, with up to 10 of the generators 
being KV. Non-conformally flat spacetimes admit a maximal Gj of conformal motions, with 
at least 1 of the generators being a proper ckv (i.e. ip ^ 0) 

Solving the conformal Killing equation (fj) to find whether, and under what conditions, a 
spacetime admits conformal motions, can lead to new solutions of the field equations (see for 
example [§] , || , H , || , ) , or to new geometrical and kinematical insights (see for example 
0' 0' @> [0)- As ^ s well-known, the standard Friedmann-Robertson- Walker (frw) 
universes are conformally flat, and admit 9 proper ckv, all pointing out of the homogeneous 
hypersurfaces M. 
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Spatially homogeneous but anisotropic universes are not conformally flat, and the exis- 
tence of proper ckv is subject to severe constraints. These are the integrability conditions 
of ([!]), including invariance of the Weyl tensor: 

A^V* = (2) 

Some conditions arise directly from ([!]) via kinematical considerations. Let w M be the four- 
velocity field orthogonal to the homogeneous hypersurfaces, so that its acceleration and 
vorticity are zero: u M = = u^ v . Then 

i> = W (3) 

C iV ? = -h^d u (u a C) (4) 

^ + (vnJ (5) 



= a 



where h^ u = g^ u + u^u u and is the shear. (The first two equations hold in any spacetime 
without acceleration and vorticity; the last holds in all spacetimes.) 

It follows from (Q) that there are no proper ckv tangent to the homogeneous hypersur- 
faces, while (|5|) shows that if is parallel to u 11 , then the shear vanishes and the spacetime 
degenerates to frw. Thus any proper ckv must be 'tilted', neither parallel nor orthogonal 
to . It is natural to impose the condition that the conformal motion of the metric should 
also be a conformal mapping of the four-velocity field w M , i.e. that £ M should map li^-curves 
into themselves [0: 

= -4)u» (6) 

However, as shown in ||, there are no perfect fluid spatially homogeneous anisotropic space- 
times with a proper ckv that satisfies the symmetry inheritance condition @. 

Thus we are led to search for a proper ckv that is tilted and does not satisfy the simple 
inheriting condition @. The simplest geometrical generalisation of (|6|) is that £ M be surface- 
forming with u^. In section 2, we impose this requirement, together with the requirement 
that £ M forms a Lie algebra with the KV of the spacetime. Given the complexity of the 
problem for general spatially homogeneous spacetimes, we begin by tackling the simplest 
anisotropic generalisation of frw spacetimes, the Bianchi I spacetimes. Then the geometrical 
and kinematical requirements allow us to solve the conformal Killing equation (|I|). The 
integrability conditions force the metric to be axisymmetric, and impose a second-order 
equation on the metric components. 

For a perfect fluid, this means that there is no freedom to specify an equation of state, 
and in section 3 we show that all perfect fluid solutions which satisfy the weak energy 
condition (p > 0) violate the dominant energy condition (p+p > 0). We consider fluids with 
anisotropic pressure, and find new radiation solutions that isotropise at late times. This 
represents a possible early-universe model with conformal motion. 

2. Solutions of the conformal Killing equation 

The diagonal Bianchi I metrics have the form 

ds 2 = -<H? + Ai{t) 2 (doff (7) 
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where x M = (x°,x l ) = (t,x), and the canonical four-velocity is = <5 M o- They admit a G% 
of motions, generated by the KV 

Y i =Y i "d li = d i with \Y i ,Y j ]=0 (8) 

From the discussion in Section 1, we are led to look for a ckv £ M that (a) closes with {Yj} 
to generate a G4, and (b) is surface-forming with u^: 

[Yi, £] = «,£ ■ 1,,'Y, (9) 
= -^i"-A/iV (10) 

where dj, are constants, and A is some scalar field, by Frobenius' theorem and (§|). 
From (P|) and (|9|) we get 

^° = a^° and = o^' + h j (11) 

The integrability condition of the second of equations (|TTJ) is 

d[jd k ]C = =>• ct[fc&j] 1 = =>- 6/ = a,£> 1 



for some constant b % (square brackets denote anti-symmetrisation). Thus (|il|) integrate to 
give 

~~'B(t)do + C\t)d t (12) 



e a-x 



where we have used the freedom to add multiples of the Yj to £ in order to set b l = 0. 
Then (0) and flT0|) imply 

C l {t) = C(t)c l (13) 
if) = Be ss (14) 

A = § (15) 

where c l are constants. To avoid degenerate cases, we require BC 7^ 0. Finally, the geo- 
metrically and kinematically defined CKV given by (0) - (|l|) is subject to the conformal 
Killing equation ([!]) for the metric (0). It is straightforward to show that the solution and 
integrability conditions are: 

axB-c l (Ax) 2 C = (16) 

aiC — = — 17 

B B A x y ! 

A 2 = A 3 = B (18) 
a 2 = a 3 = c 2 = c 3 = (19) 

It follows that such a ckv can only occur in the axisymmetric class of Bianchi space- 
times, and that eliminating C from ( |T6"D and flI7|) gives a condition involving only the two 
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independent metric functions. Writing A = A\ and a = a±, we can summarise the results as: 

ds 2 = -dt 2 + A 2 dx 2 + B 2 (dy 2 + dz 2 ) (20) 



A 2 



-rft 2 + A 2 dx 2 + 5 2 (rf?/ 2 + dz 2 ) 
e ax ^Bdt + a^J ^d?jd x 
Be ax 

BAA 2 
B~A + A2 



AB 



AB 



(21) 
(22) 
(23) 



The integrability condition (^3|) satisfies the conditions for 'decomposability' [T2|], and 
thus may be reduced to 

n 



A 
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This can be integrated upon defining a new time parameter: 

dt 
A 

B ( be aT + ce~ aT ) A 



T 



(24) 
(25) 



where b, c are constants. Clearly a ^ to avoid the isotropic frw case (B oc A). 



3. Solutions of Einstein's field equations 

The class of metrics given by (|20|) , ([24]) and fl25|) is now subjected to Einstein's field equations 

- \R5» V = pu»u u + V W V + IP, (26) 

where p is the total energy density, p is the isotropic pressure, and IL^ is the trace-free 
anisotropic pressure tensor. By symmetry there is no heat flux, and the anisotropic pressure 
has the form 

n^ = dia g (o,n,-in,-|n) (27) 

For a perfect fluid, II = 0, but II 7^ in general, e.g. for a perfect fluid with magnetic field 
ITBl). The distribution function for a kinetic gas may be expressed as [p] 
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f(x», f) = F{t, E) + F^t, E)e» + F^t, E)e»e» + ■■■ 

where F^... are the covariant multipoles that determine the anisotropy of the distribution, 
and p M = EvtF 1 + \JE 2 — m 2 e M , so that E is the particle energy and e M is a unit vector along 
the particle 3-momentum. The covariant dipole F^ determines the heat flux (which in this 
case is zero), and the covariant quadrupole F^ v determines the anisotropic pressure fl4]| : 



m 2 ) 3/2 F flu (t,E)dE 



Thus II arises from a quadrupole anisotropy of the distribution function, and by 

2A 3/2 



n(f) 



87T 
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poo , 

/ £2 - 



rn 



T{t,E)dE where F* v = diag (0, F, -§r, -±r 
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The distribution cannot be isotropic unless the shear vanishes []T5 |, so in general II 7^ (see 
16| , [H] , ||18|| , [ |i~9|j for examples). For collision-free and relaxational (bgk) distributions, 



the shear vanishes if the dipole, quadrupole and octopole are zero, or if any 4 consecutive 
multipoles are zero, or if there are only a finite number of non-zero multipoles |20f . 

Whatever the physical source of anisotropic pressure is, the field equations (|26|) with (|2T 
and (120) are of the form 



B 2 AB 
B^ + 2 AB 
_ B_B 2 _ 
~ 2 B ~ B 2 



A 
A 



B 
B 



AB 
AB 



P 

P-§n 



(28) 
(29) 
(30) 



The system (|2 



0) is closed by (|24]), fl25|) and an equation of state when II 7^ 0. 



3.1 Perfect fluid solutions 



With II = 0, there is no freedom to specify an equation of state. Subtracting 
and using (E3) and (P^), we get 



A 2 
B 2 



1 



b 2 e 2a- 

'be aT + ce" ar 
be aT - ce~ aT 



where k is a constant. Then ([[[]), (^) and @ give 

4a 2 6c /6 2 e 2ar + c 2 e~ 2aT + be' 



P 



P 



k \ b 2 e 2aT - c 2 e~ 2aT ) 
Aa 2 bc fb 2 e 2ar + c 2 e~ 2aT + 3bc" 



k 



h 2 e 2ar 



q2 g— 2ar 



from (|30| ) 

(31) 
(32) 

(33) 
(34) 



It follows from (^) and ([34]) that if p > 0, then p + p < 0, so that the dominant energy 
condition is violated, and the energy density grows with expansion, by the energy conserva- 
tion equation p + (p + p)w M ;M = 0. The limiting case p + p = occurs if and only if b or c 
vanishes, which in turn implies p = = p by (0), Q34]) , and ^ = by (|32|), fl2~2"|). In this 
case £ M degenerates to a null KV of a vacuum Bianchi I solution. Taking c = 0, we have 



ds 2 



k 

-dt 2 + a 2 (t - tA 2 dx 2 + -n:(dy 2 + dz 2 ) 

b 1 



Vk 



d t + 



1 



a(t„ — t) 



9, 



(35) 
(36) 



where t* is constant. 

In summary: although perfect fluid solutions exist for CKV that are surface-forming 
with u 11 , and do not exist at all for 'inheriting' CKV (§) that map conformally |§, these 
solutions are unphysical. 
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3.2 Anisotropic pressure solutions 

When II ^ 0, we impose the linear barotropic equation of state on the isotropic pressure: 

P = (7 " 1)P (37) 

where the constant 7 satisfies 1 < 7 < ~ for 'normal' fluids. By (^) - (|30D and ([37]), and 
using fl2"4|), (p5|), we get an equation for A(t): 



A" 



A 



12 



+ K 3 7 -4)— + 2a(3 7 -l 



+V(3 7 -2) 



A 2 

6e ar - ce~ aT 
be aT + ce _ar 



be aT - ce~ aT 
be aT + ce~ aT 



A[_ 

T 



+ 2a' 



(38) 



Thus we have reduced the problem to solving a single second order (but non-autonomous) 
ODE. We will not consider this equation in general, but only treat the autonomous case 
c = 0, when fl38D becomes 



3 4" + i (37 " 4 % + 2a(37 _ + ^ a2(37 + 2 ) = 

Although the exact solution for 1 < 7 < f has been found pill, we will only consider the 
more physically interesting case of 7 = 3. This is a radiation universe with anisotropy in 
the distribution of radiation, and ( p9|) and (E5), (|25|) give 



A 2 
T 2 



(39) 



A = fc(r — r*)e 



5 = fc6(r - r*; 



and 



a" 1 



^ e _ e -«r {a(r + 
where k,t*,r* are constants. Then (pEOD with (|2"8"D, ( |Z9"D implies 

4e 2aT [3 - 2a(r - r*)] 



n 



3fc 2 (r - 
10a(r - r«) + 3" 
2a(r — r*) — 3 _ 



P_ 
12 



(40) 
(41) 

(42) 
(43) 



4. Concluding remarks 

The exact solution (|40|) - ( f43f ) is an axisymmetric Bianchi I radiation universe with anisotropy 
in the radiation distribution, that admits a surface-forming conformal motion. There is a 
big bang at r = r#, equivalently t = t*, and for a < the universe expands as r and t 
increase. The energy density and the anisotropic pressure decay with expansion. The shear 
is given by 

1 



A^_&_ 
~A~~B 



oc 



T 



T* 



(44) 
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and decays with expansion, so that the solution isotropises at late times. The anisotropic 
pressure starts from — oo at the big bang, increases through zero, at time r = t* — (3/10a), 
to a maximum, and then decays like at late times. 

Thus the solution has some reasonable physical properties so that it might in principle 
be used as a model of the early universe. However, we are not putting it forward as such a 
model. Rather, we have shown how the imposition of natural geometrical and kinematical 
assumptions on a conformal motion lead to solutions in Bianchi I spacetimes, including 
solutions which are not immediately ruled out physically - as is often the case with conformal 
motions. 

The methods introduced here could be applied to other spacetimes, including other 
Bianchi types. 
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